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1 Introduction 



There is a deep relationship between completely integrable multi-dimensional quantum systems 
and one-dimensional multi-parameter spectral equations. One way to explain this relationship 
is to start with completely integrable models arising in the famous r-matrix approach (see e.g. 
0, III, m and references therein). As it was demonstrated by Sklyanin, many of these models 
allow a complete separation of variables in some generalized coordinate systems which, generally, 
cannot be built by means of purely coordinate transformations (see e.g. and references 
therein). The result of the separation of variables is, as a rule, a system of one-dimensional 
differential (or pseudo-differential) equations with several spectral parameters. One of these 
parameters is the energy of the initial system, while the others play the role of so-called separation 
constants. The spectral problem for the initial model reduces then to the problem of finding 
the admissible values for these parameters for which the separated equations have solutions of a 
special form. This form is dictated by the structure of the Hilbert space of the initial model and 
can be considered as an ansatz for the resulting equations. The general structure of this ansatz 
is closely related to the famous Bethe Ansatz which can be used for constructing solutions of the 
initial completely integrable system. It also turns out that the numerical equations determining 
the spectra of the multi-parameter spectral equations exactly coincide with the so-called Bethe 
Ansatz equations determining the spectrum of the initial model. 

These facts suggest to formulate an independent approach for building completely integrable 
and exactly solvable quantum systems. The idea is very simple: One should start with a certain 
A^-parameter spectral equation having an exactly constructable set of solutions and interpret N 
copies of it (rewritten in A^ different variables) as the result of separation of variables in a certain 
A^-dimensional quantum model. Identifying one of the A^ spectral parameters with the energy 
of this A^-dimensional model and considering the remaining A^ — 1 parameters as separation 
constants, one can easily reconstruct the form of its hamiltonian. Since the total number of 
spectral parameters is A^, there are A^ different ways for identifying a spectral parameter with the 
energy. This leads us to A^ independent hamiltonians whose spectra are completely determined 
by the spectrum of the initial A^-parameter spectral equation. It can also be shown that the 
resulting multi-dimensional hamiltonians commute with each other and thus form a certain 
A^-dimensional completely integrable (and exactly solvable) system. It is worth stressing that 
this (inverse) method is far from being new: It appeared approximately at the same time as 
Sklyanins (direct) method (see ref. ^) and was exposed in refs. [0, 12, 13, 14 1 where it was 
called the inverse method of separation of variables. 

With this paper we start a series of publications devoted to a detailed exposition of this 
inverse method of separation of variables. We intend to consider this method as a possible 
alternative (or rather complement) to the existing r-matrix approach. In our opinion, this 
method has several essential advantages in comparison with the r-matrix approach: First of 
all, the problem of constructing exactly solvable multi-parameter spectral equations is rather 
simple and leads to a very rich class of such equations. Applying to these equations the inverse 
method we can construct the corresponding completely integrable systems. It is easily seen 
that the class of systems obtained in this way is considerably richer than that obtained in the 
fraework of the standard r-matrix approach. The second advantage lies in the fact that we can 
consider not only exactly but also the so-called quasi-exactly solvable multi-parameter spectral 
equations. These equations are distinguished by the fact that they can be solved by purely 
algebraic methods only for some limited parts of the spectrum but not for the whole spectrum. 
The completely integrable systems obtained this way also turn out to be quasi-exactly solvable 
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(about quasi-exactly solvable models and methods for their construction see e.g. the review 
articles ||, |, |, |, [l| and the book 11] ). As far as we know, at present there are no methods 
for building such systems within the framework of the r-matrix approach. The third advantage 
finally is that the method for building exactly solvable multi-parameter spectral equations gives 
automatically the form of their solutions. This means that the problem of solving the resulting 
multi-dimensional completely integrable systems does not appear at all. We know these solutions 
from the very beginning! This is not the case for the r-matrix approach where the construction 
of explicit solutions for completely integrable models requires considerable efforts and is not 
formalised yet. 

Of course, there are also some difficulties with the inverse method of separation of variables, 
the most obvious of which can be formulated as follows: The completely integrable models 
obtained by the method are formulated in those coordinate systems in which they are separable. 
Often it is desirable to rewrite them in a form which is typical for the r-matrix approach, i.e. 
as elements of the universal enveloping algebras of some Lie algebras (we think here about 
various spin chains and their generalizations). Unfortunately, this problem is not so easy, and, 
at present, it is solved only for a few cases of models. 

2 Multi-parameter spectral equations 
2.1 General definitions 

The multi-parameter spectral equations play a fundamental role in our approach. We start this 
section giving their precise mathematical definition. 

Let Vi, i = 1, . . . , N he complex linear vector spaces of dimension di which we call carrier 
spaces. For linear operators 

eEnd[V-], i = l,...,N (2.1) 

and 

B{eEndm, j = l,...,N, i = l,...,N (2.2) 

consider the problem of finding all possible sets of complex numbers Ai, . . . , Ajv for which the 
relations 

N 

hold for some Vi 3 (pi ^ 0. We consider these relations as a system of simultaneous equations 
for the numbers Ai, . . . , Aat and vectors (pi, ... , (pjy. 

Definition 1.1. We call the system ( |2.3| ) a system of multi-parameter spectral equations. 
The numbers Ai, . . . , Aat will be called spectral parameters. 

It is convenient to regard the parameters Ai, . . . , Atv as coordinates of a point in a certain 
affine space A of dimension N. 

Definition 1.2. We call A the spectral space and those of its points for which the system 
( |2.3| ) becomes solvable we call the spectral points. The set of all spectral points will be called 
the spectrum of system (|2.3| ) and is denoted by S. 



(2.3) 
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Obviously, the spectrum T, can be found from the system of secular equations 



det,- 



A, 



N 



0, 



l,...,iV 



(2.4) 



in which detj denotes the deteminant of operators acting in the space Vi. Since (p.4|) is a system 



of equations for unknowns, the spectrum of system (2.3) is generally discrete, except for 
some special degenerate cases when it may be continuous or empty. 

Definition 1.3. We call a system of multi-parameter spectral equations finite - dimen- 
sional if the corresponding carrier spaces are all finite - dimensional. Otherwise it is infinite - 
dimensional. 

In the general case, the number of solutions of a finite-dimensional system of multi-parameter 
spectral equations is 



N. 



sol 



N 



(2.5) 



Indeed, if all the carrier spaces Vi are finite-dimensional, the secular equations ( |2.4| ) for the 
numbers Ai, . . . , Aat become algebraic equations of orders dimVi, i = 1, . . . ,N. The number of 
solutions for a system of such equations is, obviously, given by formula (^]^). If the dimensions 
of carrier spaces tend to infinity, the number of solutions of system (|2.4| ) increases. In the limit 
dimT^ = oo the spectrum of the multi-parameter spectral equations (2^) becomes infinite. This 
leads us to the following simple result: 

Proposition 1.1. Systems of finite-dimensional multi-parameter spectral equations have 
finite spectrum. The spectra of infinite-dimensional systems of multi-parameter spectral equa- 
tions are in general infinite. 



2.2 Two trivial cases 

In particular cases system ( |2.3D reduces to well known equations of linear algebra. Consider the 
following two simple examples. 

Example 1.1. Let > 1 and dimV^ = 1, i = 1, . . . , A^. Then the operators Ai and 
become numbers. The role of the vectors (pi is played by non-zero numbers which can be chosen 
arbitrarily because of the linearity of ( |2.3[) . In this case the spectral parameters Aj satisfy a 
system of inhomogeneous linear equations of the form 

TV 

A-Y,B{Xj = 0, i = l,...,N. (2.6) 

If the matrix i?^ is non-degenerate, then ( |2.6D has an unique solution 

N 

Ai = E(^"')^-4i' j = h---,N (2.7) 

i=l 

and the spectrum of system ( ^ ) consists of a single spectral point. 
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Example 1.2. Let = 1 and dimFi > 1. We denote Vi = V,Ai = A, Bl = B and also 
= (p and Ai = A. Then system ( p. 3D reduces to a single one-parameter spectral equation of 
the form 



[A-BX](j) = 0, VBcp^O. 
If the operator B is non-degenerate, then 

A G Spec \b-^A . 



(2. 



(2.9) 



In this case the properties of the spectrum of system (2.3) can be analysed simply by using the 
standard spectral theory of linear operators in vector spaces. 

The two examples demonstrate that multi-parameter spectral equations of the form ( p.3| ) 
can be regarded as generalizations of both the inhomogeneous (non-spectral) and homogeneous 
(spectral) linear equations. 



2.3 More complicated examples 

It is perhaps instructive to consider two simple examples in which the features of both the 
inhomogeneous and homogeneous equations are simultaneously present. Let us first consider an 
example of finite-dimensional multi-parameter spectral equations. 

Example 1.3. Let N = 2 and dimVi = dimV2 = 2. In this case all the operators 
Ai, Bl, B\ and B\, B2 can be represented by 2 x 2 matrices acting on two-dimensional 
vectors (pi and <j)2- The two spectral parameters Ai and A2 can be found from the system of two 
secular equations 



det (^1 - BlXi 



0, det (^2 - -B^Ai 



-B|A2 



0, 



(2.10) 



which both are algebraic equations of order two. Hence, system (2.10) has four solutions, and 
therefore the spectrum of the corresponding multi-parameter spectral equations consists of four 
spectral points. The construction of the vectors </>i and </>2 associated with these spectral points 
can be performed in a standard way. 

The next example concerns an infinite-dimensional system of multi-parameter spectral equa- 
tions. 

Example 1.4. Let N = 2. Let Vi and V2 be the spaces of analytic functions of a real 
variable x vanishing at the ends of the intervals [ai, bi] and [02, 62], respectively. Then, obviously, 
dimVi = dimV2 = 00. Let us take 



Ai = Ao 



B', 



ujI{x) 



hJ 



1,2 



(2.11) 



where loI (x) are certain analytic functions of x on the corresponding intervals. In this case, the 
system (^) takes the form 



dx'' 



w|(x)Ai - (x)A2 4>i{x)=0, 

UjI{x)Xi - UJ2ix)h <p2{x) = 0, 



h{ai) 
t>2(a2) 



Mh) 

02 (&2) 



0, 
0. 



(2.12) 
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The spectrum of this system can be found from the system of numerical equations 



^i(Ai, A2,ci,ai) = 0, 
$2(Ai,A2, 02,02) = 0, 



$i(Ai,A2,ci,6i) = 0, 
$2(Ai,A2, 02,62) = 0, 



(2.13) 



in which <I>i(Ai, A2, ci,x) and ^2(^1, A2, C2, 2;) denote the general (normalized) solutions of equa- 
tions ( 2.12| ). The system ( p.l3[ ) is transcendental and generally has an infinite and discrete set 
of solutions. Thus the spectrum of system ( 2.12 ) is infinite and discrete. 



2.4 Equivalence transformations 

Let Si E End [Vi] and Ri G End [Vi] be arbitrary invertible operators, be the elements of 
an arbitrary non-degenerate N x N complex matrix and /ij be arbitrary complex numbers. It 
is not difficult to see that the transformations of operators Ai and i?^ defined as 



Ai — Si 



N 



k=l 



RT^ G End [V^i 



l,...,iV, 



and 



= Si 



■ N 
.k=l 



R-^ G End [Vi\, 



j = l,...,N, i = l,...,N, 



and the simultaneous transformations of the spectral parameters Aj 



N 



k=l 



and vectors 



Ri 



leave the form of equations (|2.3D unchanged: 



N 



0, Vi^(t>^y^ 0, 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



Definition 1.4. The transformations (|2.14|) ~ (2.17) will be called equivalence transforma- 
tions of system (|2.3|) . The systems ( p.3[ ) and (|2.18|) related by such transformations will be called 
equivalent. We shall distinguish three particular cases of such transformations: 

1) transformations of first kind: Si 7^ li, Ri = li, Mf = 5f, /ij = 0, i = 1, . . . 

2) transformations of second kind: Si = li, Ri 7^ li, Mf = 6^, fii = 0, i = I, 

3) transformations of third kind: Si = li, Ri = li, Mf ^ 5^ , fii ^ 0, i = 1, . . 
where li denotes the unit operator acting in the carrier space Vi. 



..,N, 
,N. 



It can be seen that the transformations of the first kind do not change the solutions of 
system ( p. 3D , while the transformations of second and third kind transform the sets of spectral 
parameters Xj and vectors (pi, respectively. Note also that the transformations of second kind 
are ordinary linear transformations in the carrier spaces Vi, while the transformations of third 
kind are affine transformations in the spectral space A. 
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2.5 Degeneracy of the spectrum S 

As noted above, the spectrum of multi-parameter spectral equations can be regarded as a set of 
points in a A^-dimensional affine space A. The form of this set strongly depends on the concrete 
form of the system under consideration. In particular, it may happen that the spectrum of a 
certain system can be imbedded into another affine space A' of smaller dimension N' < N. In 
this case we shall say that this system of multi-parameter spectral equations has a completely 
degenerate spectrum. We can also imagine the situation where not the whole spectrum of a 
system but only some part of it can be imbedded into a space A' of smaller dimension. Such 
a spectrum will be called partially degenerate. Below we give a little more precise definition of 
degeneracy. 

Definition 1.5. Let S' C S denote some subset of the spectrum of some multi-parameter 
spectral equations. We call this subset degenerate if there exists a rectangular K x N matrix 

, i = N — K + 1, . . . , N, j = 1, . . . , N with K < N linearly independent rows and a K- 
dimensional vector ^j, i = N — K + 1,...,N, such that 

N 

^M/(Aj-/ii) = 0, i = N - K + 1,...,N (2.19) 
i=i 

for any (Ai, . . . , Aat) G S'. The number K (which is equal to the rank of the matrix M/) will be 
called the degree of the degeneracy. We call a degenerate subset S' non- extendable if there are 
no other elements of the spectrum S which satisfy relations ( p. 19 ) with the same M- and /ij. 



It is not difficult to see that formula (2.19) describes part of an equivalence transformation 



of third kind described in the previous subsection (see formula ( 2.16| )). This enables one to see 



that the degeneracy of a set S' C S implies the existence of equivalence transformations of third 
kind which bring the values of K spectral parameters on E' C S to zero. 

Definition 1.6. We call the spectrum S of a system of multi-parameter spectral equations 
completely degenerate if S itself is degenerate and partially degenerate if S contains a finite, 
degenerate and non-extendable subset S'. 

2.6 Hermitian multi-parameter spectral equations 

Assume now that the carrier spaces Vi are Hilbert spaces endowed with scalar products ( , )i. 

Definition 1.8. We call a system of multi-parameter spectral equations hermitian if, for all 
i = 1, . . . , A^, the operators Ai and are hermitian in the corresponding spaces Vi. 

Among many nice properties of systems of hermitian multi-parameter spectral equations 
there are two which play an especially important role in many applications. In order to derive 
them assume system ( p.3| ) to be hermitian. Then, taking the scalar product of the ith equation 
of this system with (pi, we obtain a system of ordinary linear inhomogeneous equations 

N 

Ai-J2BiXi = 0, i = l,...,N (2.20) 

with Ai = {(pi, Ai(l)i)i and = ((pi, Bj(l)i)i. Hermiticity of the operators Ai and Bj implies the 
numbers Ai and B^ to be real. But this means that the solutions Xj of system ( 2.20| ) also must 
be real. Thus we arrive at the following result: 
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Proposition 1.3. The spectra of the systems of hermitian multi-parameter spectral equa- 
tions are real. 

Let (j)f''^ and 4>(^\ i = 1, . . . ,N he two different solutions of system ( |2.3| ) with the corre- 
sponding sets of spectral parameters A^"^ and A^'"\ j = 1, . . . ,N . Then one has 



N 



i = 1, 



(2.21) 



and 



N 



0, C'^eV-, i = l,...,N 
Taking the scalar product of ( 2.21| ) with (p^-^^ and ( 2.22| ) with </)j"'', we obtain 



(2.22) 



(2.23) 



and 



m) 



(2.24) 



Due to hermiticity of the operators Ai and in the spaces Vi, we get by subtracting ( |2.23| ) 
from (|231 ) 



TV 



l,...,iV. 



(2.25) 



Let us now assume that the spectral parameters are not degenerate. Then equations ( 2.25| ) can 
be satisfied if and only if 



det 



N 



0. 



(n) 

Definition 1.9. We shall say that two solutions <j)\ and 



Am) 



(2.26) 



1, . . . , of system 



( |2.3| ) are orthogonal (in the generalized sense) if they satisfy condition ( p. 26 ) which we call the 
generalized orthogonality condition. 

Proposition 1.4. The solutions of a system of hermitian multi-parameter spectral equations 
corresponding to different spectral points are orthogonal in the generalized sense. 

It is easily seen that Propositions 1.3 and 1.4 are natural generalizations of well known 
results in linear algebra about reality of spectra of hermitian operators in Hilbert spaces and 
orthogonality of eigensolutions corresponding to different eigenvalues. 
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3 The inverse method of separation of variables 
3.1 Preliminary steps 

We will show that any multi-parameter spectral equations can be reduced to a system of ordinary 
spectral equations with a single spectral parameter. This can be done by means of the so-called 
inverse method of separation of variables {see e.g. refs. []). To explain this consider the composite 
carrier space 



V = 



Vn 



and introduce the operators 
respectively 



• Ai(S}...(S)lN GEnd [V], 



l,...,N 



j = l,...,N, i = l, 



. O/tv G End [V], 

acting in V. Here as before Ij G End [Vj] denotes the unit operator in V^. Taking 

(/? = 01 (g) . . . (g) (/)iV G V, 
we can rewrite system (^^) in the form 



N 



N 



ip = 0, V3ip^0, 



1,...,N. 



(3.1) 
(3.2) 
(3.3) 
(3.4) 

(3.5) 



This means that any solution of the initial system (|2.3D generates a solution of system (3.5). 

Now consider the N x N matrix Bj in equation ( |3.5[ ). From its definition (|3.3| ) it follows 
that the entries of this matrix belonging to different rows labeled by the index i commute with 
each other. Hence, the matrix determinant of this operator-valued matrix 



N 

P= E ^ 

jiv,ijv=i 



(3.6) 



is well defined and does not depend on the order of the factors in formula ( |3.6D . For this 
determinant, which defines obviously an operator, one can write 



N 



6ip = Y.mBi 



with 



AS]. 



TV 

E 

iivj'i-iiii+ivijjv=i 



i=l 



(3.7) 



(3.^ 



the cofactors of the elements B^. Acting by this operator matrix ( |3.8[ ) on both sides of (^^), 
and introducing the notations 



N 



1, 



.N 



(3.9) 



1=1 
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we obtain 



Cj(p = Xjpip, 9? G V, 



(3.10) 



Hence, we arrive at the following result. 

Proposition 2.1. Any solution of the system of multi-parameter spectral equations 
generates a certain solution of the system of one-parameter spectral equations ( 3.10| ). 

It is not difficult to see that the transition from ( p.SD to ( [j.lO ) essentially solves the in- 
verse problem of separation of variables in a very abstract formulation. The system of "one- 
dimensional" A^-parameter spectral equations (|2.3D is interpreted as a system of equations ap- 
pearing after the separation of "variables" in a certain "A'^-dimensional" one-parameter spectral 



equation. One of the A*" spectral parameters entering system ( |2.3| ), say parameter Aj, is identified 
with a single spectral parameter of this equation. The A^ — 1 remaining spectral parameters 
j take the role of separation constants. The reconstruction of the form of the separable 
equation is performed by eliminating the separation constants from system ( p. 31) . Since the 
index j may take A^ different values from 1 to A^, we, in fact, can associate with system ( |2.3| ) A^ 
different separable equations which are listed in formula ( ^.lOl) . 

Definition 2.1. The equations ( 3.10| ) will be called the universal separable equations. The 
operator p in this equation will be called the weight operator. 

The universality of equations ( p. 101 ) is reflected by the fact that they can be derived for 
arbitrary systems of multi-parameter spectral equations ( p.3| ) irrespective of the concrete form 
and properties of the operators Ai and Bj forming this system. 

3.2 Transformation properties of the universal separable equations 



In this subsection we derive the form of the universal separable equations ( 3.1C| ) constructed 
from the transformed equations ( ^.18| ). Introducing the operators 

5 = 5iO...(g)57v GEnd [V], 7^ = i^l®...®i^^ GEnd [V], (3.11) 

and using formulas (|3.6D - (|3.9|) , it is not difficult to check that the transformed version of system 



(|3.10D reads 



Cj(p = Xjpip, G V, j = 1, . . . , A^ 



where 



Cj=S 



N 



i,k=l 



and 



respectively 



p = 5p7^"^ 

N 



k=l 



Here we denoted by I the unit operator in the composite carrier space V. 



(3.12) 

(3.13) 
(3.14) 

(3.15) 
(3.16) 
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3.3 Transition to completely integrable models 

Up to now we did not assume any special properties for the operators Ai and Bj in ( |2.1[) and 
( |2.2| ). In this section we show, that there exists a wide class of these operators for which the 
corresponding universal separable equations (|3.10| ) can be interpreted as spectral equations for 
a certain completely integrable quantum system. We start with the following definition: 

Definition 2.2. We call the system of multi - parameter spectral equations non - degenerate 
if its weight operator p defined by formula (3.6) is non-singular, i.e. if det p 7^ 0. 



Invertibility of the weight operator implies the existence of operators 

Vi = p-'^d £ End [V], i = l,...,r (3.17) 

whose eigenvalues in the space V coincide with the values of the parameters Aj: 

V,^ = Xj^, 99 eV, j = l,...,N (3.18) 

From ( |3.18| ) it immediately follows that 

[Vj,Vk]ip = 0, ipeV, j,k = l,...,N (3.19) 

for any ip satisfying equations (|3.1 j ). This means that the operators Vj commute with each 
other in the week sence, i.e. on the states ip built from the solutions of the initial system (p.3|). 
Commutativity in the strong sense does not follow generally from ( 3.19| ) if it is not a priori 
known that the set of these states is complete in V. 

Fortunately, there exists a special case in which strong commutativity of operators Vj can 
be guaranteed without any knowledge of the completeness of ip^s. This case happens when the 
operators Bj, j = 1, . . . , N commute with each other for fixed i = 1, . . . , N: 

[Bj,B^]=0, i,j,k = l,...,N. (3.20) 

Definition 2.3. We call a system of multi-parameter spectral equations commutative if the 
operators B^- in (|2.2|) satisfy relations ( [3.20 ). 



In order to demonstrate that commutativity of system ( |2.3D implies commutativity of the 
operators Vi in ( ^.17 ), let us first note that if the operator p is invertible, the operators Vj can 



be defined as solutions of the system of linear operator equations 



N 

Ai = Y.B{Vj, i = l,...,N. (3.21) 



From definitions ( |3.2D and (|3.3| ) of the operators Ai and and formula ( 3.20 ) it follows that 



[Ai,Ak]=0, for any i,k, (3.22) 



0, for all i, k,j, I 



(3.23) 



[A. 



0, for any i ^ k and any j. 



(3.24) 
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Taking into account relation ( 3.2l| ), we get the following chain of equalities: 



N N N 

n,m=l n,m=l n,m=l 

N N N 

+ E [B^,BT]VnV^+ E B^BT[Vn,Vj^= E [B':Vn,BTVm] = 

n,m=l n,m=l n,m=l 

= [Ai,Ai]=0. (3.25) 

Here we used the fact that the last two terms in the first line of ( p. 251) c ancel and the first term 
in the second line vanishes because of the commutation relation ( p.23|) . Now, using relations 
(3.22) and ( 3.2^ ), we can consider an analogous chain for i ^ k: 



N N N 

E B?BnT^n,V^]= E B?Bn'DnM+Y.i'3'^l3T]'^m = 

n,m=l n,m=l m=l 

N N N 

E B?Br[Vn,nm]+ E [B'^'Dn,Br]v^= E B?Bn'Dn,Vm] + 

n,m=l n,m=l n,m=l 

N N N 

+ E B^[Vn,B^]Vm= E Bnv^,BtV^] = Y.m^n,Au] + 

n,m=l n,m=l n=l 

N N 



+ E [K, Bl]Vn = E [^r^n, Ak] = [A, Ak] = (3.26) 



n=l n=l 

Combining ( |3.25D and ( p. 26 ), we obtain finally 



N 

E B^Bf[Vn,Vj^ = Q, for any i and fc. (3.27) 

n,m=l 

Now note that the invertibility of the operator p implies invertibility of the matrix operator 
(|3.3| ). Using this fact, we obtain the relations 

\Dn-,T>m\ = 0, for any n and m, (3.28) 

i.e. commutativity of the operators Vi in the strong sense. Thus, we proved the following 
theorem: 



Theorem 2.2. If the system of multi-parameter spectral equations ( p.3| ) is non-degenerate 
and commutative then the operators Pj, i = 1, . . . , N defined by formulas ( |3.17 ) form a com- 



mutative family and can be considered as integrals of motion of a certain completely integrable 
quantum system. These operators have a common set of eigenfunctions (/? constructed from the 
solutions of system (^]^) and their eigenvalues coincide with the values of the spectral parameters 



Aj, i = 1, . . . , of equation (|2.3|) . 



3.4 Transformation properties of completely integrable equations 

It is quite obvious that equivalence transformations do not affect the invertibility properties of 
the operator p, as follows from formula ( 3.14| ). However, these transformations may change the 
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commutativity properties of the operators Bj . In other words, the commutativity of a certain 
system of multi-parameter spectral equations does not automatically imply the commutativity 
of its transformed version. However, it turns out that the transformed version of equations 
( 3.18 ) describes again a certain completely integrable system. Indeed, applying to ( p.l8 ) some 
equivalence transformation, we obtain 



where the transformed operators Vj, using formulas ( 3.17D , (|3.13|) and ( ^.14 ), read as 



Do = 7^ 



N 



J2 M^{Vk - iMkI) 



.k=l 



7^" 



But from ( p.28| ) and ( p. 30 ) it immediately follows that for any n and m 



(3.29) 



(3.30) 



(3.31) 



i.e. the transformed model is again integrable. This leads us to an important conclusion which 
we formulate in the form of the following theorem: 

Theorem 2.3. Let a system of multi-parameter spectral equations be non-degenerate and 
equivalent to a certain commutative system. Then the operators Vi, i = 1, . . . ,N defined by 
formulas ( |3.17 ) form a commutative family and can be considered as integrals of motion of a 
certain completely integrable quantum system. 



3.5 Hermitian completely integrable systems 

Let us now consider a system of hermitian multi-parameter spectral equations. Remember that 
hermiticity means that the spaces Vi are endowed with scalar products ( , )i with respect to 
which the operators Ai and Bf are hermitian in Vi. Then formulas ( |3.lD allow one to introduce 
in a natural way a scalar product ( , ) in the space V. The operators Bj will obviously be 
hermitian in V. Using the fact that the operator Bi commutes with all the operators Bl. with 
k ^ i and observing the absence of jB" in expression ( |3.8D , we can conclude 



/^B),Bi 



■ 0, for any j and i. 



(3.32) 



Note also that, due to the commutativity of the operators Bj and their hermiticity in Vi, the 
operators (|3.6D and (|3.8|) are hermitian in V: 



(det B)+ = det B, {B)r = B]. (3.33) 
From ( |3.32| ) and ( 3.33 ) it follows that the operators Cj and p are hermitian in the space V. 



Definition 2.4. We call a system of multi-parameter spectral equations positive definite if 
it has a positive definite weight operators. 

It is obvious that positive definiteness of a hermitian operator p implies its invertibility and 
also the existence of the square root p^/^ which can also be considered an hermitian operator. 
Let /3 be a positive definite hermitian operator. Then, after introducing new vectors ^ as 



(3.34) 
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and new operators Dj as 

equation ( |3.10| ) takes the form 

Vj^ = Xj^, cI>eV, 3 = 1,..., N 



(3.35) 
(3.36) 



The operators Vj are obviously hermitian in the space V and their eigenvalues in V coincide with 
the admissible values of the spectral parameters \j. It is not difficult to see that the operators 
(3.35) are related to the operators introduced in ( 3.17] ) by the formula 



P 



(3.37) 



From this formula it follows that commutativity of the operators Vj implies commutativity of 
the operators V. But from previous results we know that commutativity of the operators T) is 
implied by commutativity of the initial system of multi-parameter spectral equations. Collecting 
all these facts and recalling the definition of exact and quasi-exact solvability we arrive at the 
following important theorem: 

Theorem 2.4. Let a system of multi-parameter spectral equations be hermitian, positive 
definite and commutative. Then the operators Vj, j = 1, . . . , N are hermitian, form a commu- 
tative family, [VjjDk] = 0, and thus, can be considered as "hamiltonians" of a certain comletely 
integrable quantum system. If, in addition, the system of multi-parameter spectral equations 
is exactly or quasi-exactly solvable, then the resulting completely integrable system also will be 
exactly or quasi-exactly solvable. 

This theorem will play a central role in our further considerations. Before completing this 
section, it is important to discuss the case when the initial multi-parameter spectral equation 
becomes degenerate. 

3.6 The degenerate case 

Assume that the operator p is non-negative definite and hence singular in V. In this case, for 
any j = 1, . . . ,N there exists an orthogonal decomposition of the space V of the form 

V 



in which the operators Cj and p take the following block form 



(3.38) 







nl2 
^3 


^3 






/r 









n2\ 


n22 
^3 







p = 













\cf 








) 




[ 





J 



(3.39) 



Here r E End[V^] is a positive definite hermitian operator, Cj^ € End[V^] and G End[V?] are 
hermitian operators, and € Hom[V|, V^] and G Hom[V|, V^] are hermitian conjugates of 
the operators Cj^ G Hom[V^, V|] and Cj"^ G Hom[V^, V|]. The decomposition ( |3.38| ) can always 
be chosen in such a way that all the operators CP are invertible. Writing 



(3.40) 
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with (f^ G , (pj Vj and € V|, we can rewrite equations ( |3.10|) in the form 






From the last equation in ( |3.41| ) it follows that 

= P^' 

where P is a certain operator such that 

CfP = 0, P^Cf = 0, 3 



(3.41) 



(3.42) 



(3.43) 



and if' € V' is a vector belonging to a certain subspace V' C V^. Substituting ( 3.42| ) into the 
second equation in (|3.41|) and using ( 3.43| ) we obtain 



^1 



{cfr^cfp^' 



(3.44) 



Substituting ( 3.42| ) and ( p.44 ) into the first equation in ( 3.41| ) and acting on the resulting relation 
by the operator P~^ we obtain the following set of one-parameter spectral equations 



/->/ / 



in which 



C' = P' 



are hermitian operators in V' and 



= Ajp if , if £ V 

'cf-cf(cf)-^cf 

p' = P+RP 



p 



(3.45) 



(3.46) 



(3.47) 



is hermitian and positive definite in V'. Noting that {p'Y^"^ is an invertible hermitian operator 
and introducing the vectors, 



and the operators 



{pr"'c'^{p')-''\ 

we can reduce equations ( |3.45| ) to the form 

Xj^', 1>'gV', j = l,...,iV 



(3.48) 



(3.49) 



(3.50) 



The operators V'^ are obviously hermitian in the space V and their eigenvalues in V' coincide 
with the admissible values of the spectral parameters \j. As before, if the set of solutions is 
complete in V', the operators T>'j, j = 1, . . . , N form a commutative family, [P^jP^ = 0, and 
thus, can be considered as "hamiltonians" of a certain completely integrable quantum system. 
Unfortunately, at present we do not have more general criteria for commutativity of the operators 
Vj in the degenerate case. 
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4 Conclusion. Exact and quasi-exact solvability 



Let us now discuss some solvability properties of systems of multi-parameter spectral equations. 
First note that if all the carrier spaces Vi are finite-dimensional, then the solution of the system 
of multi-parameter spectral equations is a purely algebraic problem (see e.g. examples 1.1 and 
1.3). In the case of infinite-dimensional carrier spaces an algebraic solution of the system of 
multi-parameter spectral equations is, as a rule, impossible (see e.g. example 1.4). However, the 
presence of some high symmetry in the system may turn it algebraically solvable, despite the 
infinite-dimensionality of the carrier spaces. 

Definition 4.1. We call an infinite-dimensional system of multi-parameter spectral equa- 
tions exactly solvable if all its solutions can be constructed algebraically. We call an infinite- 
dimensional system of multi-parameter spectral equations quasi- exactly solvable if only a certain 
finite part of its solutions allows an algebraic construction. 



Consider an infinite-dimensional system of A^-parameter spectral equations ( |2.3| ) which is 
known to be exactly solvable. Assume that the spectrum of this system is partially degenerate 
and thus, by definition, contains at least one finite, degenerate and non-extendable subset E'. 
Denoting its degree of degeneracy by one can write 



N 



0, 9 0, ^0, (Ai,...,A^) G S', i = l,...,iV. 



(4.1) 



According to the results of the previous subsection, the degeneracy allows to apply to relations 
(|4.l| ) an equivalence transformation, which turns them into the form 



N' 



0, V^3^i^Q, (Ai,...,Aiv') e S', i = l,...,iV', 



(4.2) 



Here we used Aj = for j = — -|- 1, . . . , A^, suppressed the last K equations for i = 
N — K + 1, . . . , N and introduced the notation N' = N — K. We also denoted by S' the set 
of the transformed spectral points with the zero components suppressed. It is easily seen that 
system ( |4.2| ) is nothing else but the system of A^'-parameter spectral equations whose spectrum 
contains by construction a certain finite exactly solvable part S'. Since the tranformed system 
( |4.2| ) is infinite-dimensional, its total spectrum S is infinite. This means that the A^'-parameter 
spectral equation obtained this way is quasi-exactly solvable. Thus we arrive at the following 
theorem: 

Theorem 4.2. Any infinite-dimensional and exactly solvable system of A^-parameter spec- 
tral equations with partially degenerate spectrum can be reduced to a certain quasi-exactly 
solvable A^'-parameter spectral equation. Here N' = N — K where K denotes the degree of the 
degeneracy of the spectrum of the initial equations. 
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